Stability of magnetic nanoparticles inside ferromagnetic nanotubes 
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During the last years great attention has been given to the encapsulation of magnetic nanoparti- 
cles. In this work we investigated the stability of small magnetic particles inside magnetic nanotubes. 
Multisegmented nanotubes were tested in order to optimize the stability of the particle inside the 
nanotubes. Our results evidenced that multisegmented nanotubes are more efficient to entrap the 
particles at temperatures up to hundreds of kelvins. 
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Interest in magnetic nanostructures has increased dra- 
matically over the past decade, mainly due to the great 
progress in experimental techniques and recent techno- 
logical developments which allow access to interesting 
length scales. In particular, understanding the behaviour 
of solids in confined geometries at the nanoscale is an 
important area of basic research which has a strong po- 
tential for applications^— . Within this area, carbon nan- 
otubes have stimulated intensive experimental and the- 
oretical research since their discovery in 1991—. One of 
the topics that have been considered is the encapsula- 
tion of nanoparticles for different purposes. In particu- 
lar, paramagnetic nanoneedles have been encapsulated 5 , 
preventing them from agglomeration when a magnetic 
field is not present and allowing control of the needles' 
movement. The carbon encapsulation of ferromagnetic 
nanoparticles has also been reported*—, opening an in- 
teresting field of research. In particular, encapsulation of 
NiS, CoiS, Feii~— and Ga^ inside multi-walled carbon 
nanotubes has been reported and magnetically charac- 
terized using vibrating sample magnetometers and su- 
perconducting quantum interference devices. Also, the 
stability and magnetic properties of Fe encapsulated in 
silicon nanotubes has been explored, finding a high local 
magnetic moment per Fe atom, and a transition between 
from antiferromagnetic to a ferromagnetic coupling is ob- 
tained by increasing the length of the nanotubei 5 -. 

Besides this intense study on the encapsulation of 
magnetic nanoparticles, few attention has been given 
to the encapsulation of nanoparticles inside magnetic 
nanotubes^. Magnetic nanotubes have been prepared 
since 2005^ by K. Nielsch and his group. Also barcode- 
type nanostructures have been prepared. In particular, 
Lee et a/.— reported the synthesis of multisegmented 
metallic nanotubes with a bimetallic stacking configu- 
ration along the tube axis, showing a different magnetic 
behavior as compared with continuous ones. 

Following these ideas, the focus of this work is the sta- 
bility of a magnetic nanoparticle inside a ferromagnetic 
nanotube. Using Monte Carlo simulations we investi- 
gate different nanotubes, looking for the magnetic con- 
figuration that increases the stability of an encapsulated 



magnetic nanoparticle. To investigate the behaviour of a 
nanoparticle inside a magnetic nanotube we start by cal- 
culating the magnetic fielcU^ due to the nanotube over 
the points of a three-dimensional grid ({ri}). This step 
is accomplished by adding up the contributions to the 
overall magnetic field coming from each atom that con- 
stitutes the nanotube (the magnetic dipolar field). Then 
we simulate the random walk of a magnetic nanoparti- 
cle over the field landscape generated by the nanotube. 
This is done by means of a Monte Carlo simulation in 
which the nanoparticle starts at the center of the nan- 
otube and is allowed to move over the grid and to rotate 
its magnetic moment under the influence of the tube's 
field. 

The starting point is the definition of the crystalline 
structure of the tube, i. e., we define the Bravais lat- 
tice, the lattice constant (a ) and the magnetic moment 
per atom (jj, ). The simulations shown in this paper con- 
sider BCC iron as the constituting species, which gives us 
a a = 0.2866 nm and ^ = 2.22 ^ B . An important point 
is the definition of the three-dimensional grid over which 
the field will be calculated and the energy of the parti- 
cle will be tested. For our calculations we have chosen 
a grid that is commensurate with the lattice structure 
and that considers the available computational facilities. 
According to these requirements, we have used as the 
grid step (A), the nearest multiple to 5 nm of a , hence, 
A = 17 a„ = 4.8722 nm. However, we have also tried 
smaller and larger A (A = 8 a a and 52 a , respectively) 
and no fundamental differences were observed compared 
with the results obtained using A = 17 a . 

The next step is to define the tube's overall dimen- 
sions. The tube's height, which is also commensurate 
with a , was chosen to be Ltubc = 205 A = 998.8 nm. 
The radial dimensions were set to R c = 50 nm (external 
radius) and Ri = 40 nm (internal radius). We investigate 
the behaviour of the particle using multilayer nanotubes 
composed of magnetic and non-magnetic segments, as 
illustrated in Fig. 1. The length of the non-magnetic 
segment, the spacer, has also been varied (see Fig. 2) to 
obtain the desired field landscape, but similar results are 
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obtained if one uses a non optimal choice. 
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FIG. 1. (Color online) Illustration of the different tube ge- 
ometries that we have tested, depicting the magnetic and 
non-magnetic regions. The white arrows represent what we 
call an antiparallel segments. 



the tube's axis obtained from Eq. 2 for one or two seg- 
ments while keeping the full length of the tube fixed at 
998.8 nm, as illustrated in Fig. 1. In this figure we can 
observe that for thinner spacers the antiparallel align- 
ment provides a higher magnetic field in the centre of 
the tube compared to the parallel alignment. 
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Even a nanotube with such dimensions would contain 
more than 240 million atoms, rendering extremely ex- 
pensive any attempt of direct calculations, regardless of 
the computational power available. To circumvent this 
technical issue, we first calculate the field generated over 
the grid points ({r^}) by all the magnetic moments in- 
side a ring (p ) of height A which is magnetized along its 
axis and located at the centre of the grid. We used the 
usual expression (Eq. [IJ, for the field of a collection of 
magnetic dipoles and stored our results for each point of 
the grid in a file. The latter was used as a building block 
to construct the desired tube. 

z3 /->\ 3^t Q • (r j — r a ) 
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T3 + Y&«- S (fi - r a ) , (1) 



where the indices a and i run over all the atoms belonging 
to the ring p and all the sites belonging to the grid, 
respectively. 

Finally, by summing the contributions coming from 
different rings (p n ) centered at different positions (f n = 
nAz) with different weights (w n ), we are able to con- 
struct multi-segmented tubes^ with as many parallel or 
antiparallel aligned segments as we want (just by weight- 
ing their contribution by ±1) and separated by non- 
magnetic materials (by weighting their contribution by 
0) as follows. 



-Btube (fi) 
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w n B (n - f n ) , 



(2) 



where N — (L tube — A) /2A. The result of such calcu- 
lation is depicted in Fig. 2, which shows the field along 



FIG. 2. (Color online) The magnetic field along the tube's 
axis obtained from Eq. [2] for one or two segments while keep- 
ing the length fixed. 

Above we have defined the set of grid points occu- 
pied by the tube itself and the magnetic field B ({fi}) 
that it creates. So now we are in position to investi- 
gate the stability of a nanoparticle inside the tube. We 
start by putting the particle, characterized by its posi- 
tion (f £ {f}) and magnetic moment vector (/t), in the 
center of the tube and let it move and rotate under the 
Hamiltonian 
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The magnitude of the moment p is chosen to repre- 
sent a Co nanoparticle with diameter A whose mag- 
netic moment initially points in a random direction. At 
each Monte Carlo step, a movement and a new orien- 
tation for the magnetic moment are proposed simulta- 
neously. The movement consists in walking one step 
on the grid along one of the six possible directions 
(+x, — x, +y, — y, +z, — z). The new orientation is sam- 
pled randomly from a sphere with radius \fl\. 

The energy difference caused by the new position and 
orientation is calculated and the Monte Carlo step is ac- 
cepted (or not) by means of the standard Metropolis al- 
gorithm at a given temperature T— . All the simulations 
shown on this paper were done using a total of 10 7 Monte 
Carlo steps. Then we repeat this procedure 2 x 10 3 times 
in order to get the average walk of the particle on the 
magnetic field landscape of the tube. 

We have performed the above-mentioned procedure for 
different tube geometries, magnetic configurations and 
temperatures. We have chosen the geometry in which 
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the non-magnetic spacer is d — 51a thick because thin- 
ner spacers would allow long-range magnetic interactions 
(such as RKKY) to play an undesired role by strengthen- 
ing the coupling between different segments. Our results 
are depicted in Fig. 3, which shows a measure of the aver- 
age z coordinate (along the tube's axis) of the particles as 
a function of the number of Monte Carlo steps (mcs). The 
graphs are shown in terms of £ = \J (^ 2 }/(Ltube/2). £ = 
denotes that particle is at the center of the tube, £ « 0.34 
represents one of the interfaces between the spacer and 
the magnetic section in the three-segment nanotube, the 
tips of the tube are denoted by <; = 1, and a particle 
leaving the tube is denoted by £ > 1. 

Figure 3 shows that at T = 10 K all the compound 
tubes are able to confine the particle while the simple 
tube is not. At T — 300 K none of the tubes were ca- 
pable of keeping the particle trapped, although the tube 
with two antiparallel aligned segments was the one that 
mostly delayed the leaving. However, it is at T = 100 K 
that interesting phenomena take place: the parallel align- 
ment of segments seems to be less effective in trapping 
the particle, since only the field landscape created by an- 
tiparallel aligned segments is capable of restraining the 
particle's movement. 
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FIG. 3. (Color online) Evolution of f as a function of the 
number of Monte Carlo steps (mcs) for several tube geome- 
tries (one, two or three segments), different magnetic config- 
urations (parallel or antiparallel alignment) and several tem- 
peratures (T = 10, 100, 300 K). 5 = 1 represents the tips of 
the tube. 

We also investigate the effect of the size of the nanopar- 
ticle by increasing its diameter from A to 2A. Our re- 
sults for T — 300 K are presented in Fig. 4, showing that 
larger sizes contribute to stabilize the particle inside the 
tube. 

From our results we can conclude that the use of multi- 
layer geometries can help to stabilize magnetic nanoparti- 
cles inside ferromagnetic nanotubes. In spite there is still 
no experimental realization of these proposal, it might 
be possible to synthetized magnetic nanoparticles inside 
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FIG. 4. (Color online) Evolution of £ as a function of the 
number of Monte Carlo steps (mcs) for particles with diameter 
A (top) and 2 A (bottom). 



nanotubes by coating a multi-segmented nanotube with 
an oxide layer on the inner wall, using Atomic Layer De- 
position (ALD). Subsequently, using electrodeposition, 
the nano-channel can be filled up to the middle of the 
pore channel with the non-magnetic segment. At the 
end of the electrochemical deposition a short segment 
of magnetic material can be deposited. Similar tech- 
niques have been used to produce multilayered core/shell 
nanowires by Chong et a& and multi-segmented nan- 
otubes by electrochemical deposition 18 . Before starting 
the experiment, the non- magnetic metallic segment is se- 
lective etched and the oxide layer should be also selec- 
tively removed. During the removal of the oxide layer on 
the nanotube inner wall, and using a very strong mag- 
netic field in the perpendicular direction to the nanotube 
axis, it is possible to trap and mechanically detach the 
magnetic particles from the template. By removing the 
perpendicular magnetic field the time until the particle 
exit the tube can be measured. Alternative it is possible 
to make the magnetic segments of a material with a Curie 
temperature around 100 °C and inject the magnetic par- 
ticle with a polymeric meld at a higher temperature into 
the channel. After cool down the polymeric meld, it will 
solidify. The polymeric filling can be removed using an 
organic solution and a strong magnetic field perpendicu- 
lar to the nanotube's diameter. In this way, these proce- 
dures can guide the design of experiments which can test 
our results. 

Looking at what happened at T = 100 K, it seems 
possible to use the alignment between the segments as 
a switch to open or close the particle's exit. If one con- 
structs a tube with two segments, one made of a soft 
and the other of a hard magnetic material, one could 
in principle use an external magnetic field to switch the 
direction of magnetization of the softer segment, thus 
changing from antiparallel to parallel alignment and fa- 
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cilitating the departure of the walking particle. 

If one wants to increase the temperature where 
the free/trapped behaviour occurs when using paral- 
lel/antiparallel alignment (in our case, T = 100 K for 
particles with diameter A) one can simply use walking 
particles with larger magnetic moments to increase the 
magnitude of the Zeeman energy (i? Zcoman = —jl- B) with 
respect to the thermal energy (k B T). 
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